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Prataee 

These three reports describe selected applications that the Texas Instruments Digital Signal Processing 
(DSP) Gioiqi, in ooajuncdon widi H leseardi gtoapB, has develqwd in the emeigiiig field of fbzzy logic. 

The first paper gives an overview of the fiizzy logic tiwory and sugge^ methods tS. Inqdonoitation that 

are illustrated in the next two papers. 

The second paper describes the implementation of fuzzy logic as a serial algorithm on a TMS320 DSP sys- 
tem. The algorithm is based on parts of the theory described in the third paper. 

The third paper presents a Aeoiy for a dedicated siliom-level fiizzy logic implementation. 
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What is Fuzzy Logic? 
An Overview of thie 
Latest Control {Methodology 



Hmothy A. Adcock 
Digital Signal Processing — Semiconductor Group 
lexBS IiM«Qii&ls Ineoir^rated 



Introduction 



The name "Fuzzy Logic" seems to imply an imprecise methodology that is useful only when accuracy is 
not necessary or important. That is what many people assume when they first hear about fuzzy logic — and 
understandably so. In a world increasingly manipulated by computers with their absolute "1" or "0" and 
"on" or "off concepts, a term like fuzzy logic suggests inaccuracy or imprecision. Even Webster's dictio- 
nary defines '%izzy" as: 

ftizz*y (-^ a^}. 2. not clear, distinct, or precise; blurred 

This is not true of fuzzy logic. Fuzzy logic can address complex control problems, such as robotic arm 
movement, chemical or manufacturing process control, antiskid braking systems, or automobile transmis- 
sion control with more precision and accuracy, in many cases, than traditional control techniques have. 

Fuzzy logic was invented and named by Lotfi Zadeh, a professor at the University of California at Berkeley. 
Fuzzy logic is a methodology for expressing operational laws of a system in linguistic terms instead of 
mathematical equations. Many systems are too complex to model accurately, even with complex mathe- 
matical equations, but fuzzy logic's linguistic terms provide a useful method for defining the operational 
characteristics of such a system. These lingui^c tenas Me most ^ea expressed in the form of logical 
implications, such as If -Then rules: 

Ifairjemp is WARM, then setfanjq>eed to MEDIUM. 

The terms WARM and MEDIUM are actually sets that define ranges of values known as membership func- 
tions. By choosing a range of values instead of a single discrete value to define the input variable 
''air_temp", you can control the output variable ''fon_speed" more precisely. Fuzzy logic controllers can 
often improve the performance of a control system by reducing the diance of wild functions in the output 
that may be caused by variations in the measured input variables. 

The IVadltional Approach 

To illustrate the difference between fuzzy logic and the traditional approach, here is a control problem. 
First, consider how the trad^nd— called "crisgT — craitiollw would handle it: 

Ifati-Jemp is S 70" Fahrenheit, then setfan_speed to "1000 rpm". 
Ifair_tenip is < 70' Fahrenheit, then set fan_speed to "100 rpm". 

A nonfuzzy, or "crisp," ctmtroller relies on a discrete valued decision point. For this type of system, the 
input must reach an exact value before the control system reacts in a certain way. Even small variances in 
this input value may cause the output to react drastically differently. For instance, if the temperature is 70° 
or above, the first rule will set the fianjspeed to 1000 rpm. If the temperature is below 70", die second rule 
will set the £ui_speed much lower, to 100 ipm. Figure 1 shows a diagram of tiiis oisp valued controller. 
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What would happen if the temperature were 69.5 " ? Or, more importantly, what would happen to the control 
system if the temperature were transitioning from below 70° to above 70"? The temperature might even 
fluctuate back and forth slightly above or below 70° (e.g., 69.0° to 71.0°). This would cause the control 
system to alter the fan speed wildly for changes in the iapat variable air_temp, although the temperature 
change may not be significant. 

These transition points are difficult for "crisp" control systems to handle, but they are exactly where "fuzzy 
logic" excels. 

Fuzzy Control 

Fuzzy logic is implemented in tiiree phases (see Figure 2): 

1. Fuzzification (crisp input to fuzzy set mapping). 

2. biference (fuzzy rale generation). 

3. Defuzzification (fiizzy to crisp output transformation). 

Figure 2. Fuzzy Logic Phases 
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Fuzzification 

In the first fuzzy logic phase — ^fuzzification — actual measured input values are mapped into fuzzy mem- 
bership functions. As an example, a climate-control system has been developed with fuzzy logic. 

To create a climate control system, we first developed membership functions for the input variable 
"air_temp". These membership functioiis are defined by both a range of values and a degree of member- 
ship. In fuzzy logic, it is important to distinguish not only which membership functions a variable belongs 
to, but also the relative degree to which it is a member. This gives the variable a "weighted" membership 
in a membership function. A variable can have a wei^ted member^p in several membership functions 
at the same time. The membership functions ftx "aktemp" are shown ki Figure 3i 
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Fisur» 3. A(r_T«nip tii|NJt Vnrtobtle Ntambiirship Functions 




As shown in Figure 3, fuzzy membership fiinctions span a range of values and can actually overlap. Three 
sets of membership values are defined above for the variable "air_temp". They are COOL, WARM, and 
HOT. The degree of membership is found by finding the intersection point of a distinct input value on the 
horizontal axis with the line defining one or more fuzzy membership functions. This intersection point is 
assigned a corresponding value on the vertical axis to define the relative membership in a set for an actual 
measured input value. Notice that when "air_temp" is at a particular value, it may be contained in one or 
more fuzzy sets. For instance, at 70° , "air_temp" is a member of the function HOT with a relative member- 
ship of 0. 17. It is also a member of the function WARM with a relative membership of 0.37. Unlike a crisp 
system in which a value either is or is not a member of a function, a fuzzy logic system can take action based 
not only on membership in a fuzzy set, but also on the degree to which a variable is included in a member- 
ship function. In this case, because "airjemp" at 70° is more WARM (0.37) than it is HOT (0.17), the 
controller will take that into account when defining what output a<:tion to take. 



Inference Rule Definition 



Once membership functions have been defined for input and output variables, a control rule base can be 
developed to relate the output actions of the controller to the observed inputs. This phase is known as the 
inference, or rule definition portion, of fuzzy logic. Any number of rules can be created to define the 
actions of the fuzzy controller. Some examples are shown below. 



Fuzzy Logic Rule Definition 

If air_temp is COOL, then set fan_speed to SLOW. 

If air_temp is HOT, then set fan_speed to FAST. 

If air_temp is WARM, then set fan_^ed to MEDIUM. 
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These If-Then rules can relate multiple input and outpmt variables. Because the rules are based on ■ward 
descriptions instead of mathematical definitions, any relationship that can be described with linguistic 
terms can typically be defined by a fuzzy logic controller. This means that even nonlinear systems can be 
described and easily cmtrolled with a fuzzy logic controller. In addition, since variable have weighted 
memberships — in particular membership functions — the rules that are composed of these variables are 
weighted as well. This means that different rules have different impacts on the controller, according to the 
measured input variable. For a multiple-input/multiple-output system with many defining rules, a wild 
fluctuation in any single input will be tempered by these rule weightings. Because of this, fuzzy logic sys- 
tems are very robust and often allow many rules to be removed or altered without significantly impacting 
the controller. 

Defuzzificatlon 

After the fuzzy logic controller evaluates inputs and applies them to the rule base, it must generate a usable 
output to the system it is controlling. This may mean setting a voltage or current to a particular value to 
control the speed of a fan in the example above, or it may mean defining the optimal speed of a robotic arm 
as it nears its target. The fuzzy logic controller must convert its internal fuzzy output variables into crisp 
values that can actually be used by the controlled system. You can perform this portion of the fuzzy control 
algorithm, known as defuzzificatlon, in several ways. IVo of the most common methods are: 

• maximum defuzzification method (page 6). 

* centroid calculation defuzzification method (page 7). 

Remember from fuzzification that in mapping input variables to membership functions, a particular mea- 
sured value of the input variable determined the relative membership of that input variable in an input mem- 
bership function. To determine tiie mapping of output variables to tiieir corresponding output membership 
functions, the weighted input membership function and corresponding rule base determine the relative 
membership in the output function. Whatever relative membership was given to the input variable will also 
be given to the output variable, as assigned by its corresponding rule. For air_temp = 70° , the output vari^ 
aides are assigned a value that corresponds to the input value shown in Table 1. 

Table 1. Fan_Speed (Membership Function Relative Membership) 



Input Variable 


Defining Rules 


Output VariaUe 


air_temp (WARM) = 0.37 


If airjemp = WARM, 

then set fan_speed to MEDIUM 


fan_speed (MED) = 0.37 


airjemp (HOT) = 0.17 


If airjemp = HOT, 

then set fan_speed to FAST 


fan_speed (FAST) = 0.17 



The output variable fan_speed is given the same relative mapping as the input variable air_temp that is 
defined by a particular rule. 

Figure 4 illustrates the output variable membership functions. It this case, a distinct value on the horizontal 
axis is defined by the relative membership on the vertical axis. To create tiie actual crisp output value for 
the controller system output, membership functions are used with the output variable. In Table 1, the input 
value air_temp = 70 resulted in two weightings for fan_speed: 

* Fan_speed = 0.37 was assigned to the output membership function MEDIUM. 

* Fan_speed - 0.17 was assigned to the output membership function FAST. 
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As shown in Figure 4, the actual output value is detennined by beginning at the weighting £actor on the 

vertical axis and moving horizontally until an intersection point is reached on the lines defining its 
associated membership function. This intersection point is then transposed to the horizontal axis to deter- 
mine the crisp output value. 




Fan Speed (rpm) 



Maximum Defuzzlflcatlon Method 

One method of defuzzification is known as the maximum method. In this method, if more than one rule 
is active, the maximum relative membership is used to determine die output vtioe. hi the above example, 
maldng air_temp = 70° created two possible values for fim^speed: 

fan_speed = 0.37 
fim_speed = 0.17 

The maximum defuzzification method provides a single output by choosing the active rule with the greatest 
relative membership value in the output membersMp^ction. In the preceding example, the following rule 
is chosen because it has the highest membership value for fan_speed. 

Ifairjemp = WARM, then setfanjspeed to MEDIUM. Fan_speed (MED) = 0.37 

The value 0.37 on the vertical axis intersects the membership function MEDIUM at two points — one on 
the positive slope (at 375 rpm) and one on the negative slope (at 710 rpm) of the function. The two points 
represent two possible solutions that must be resolved. 
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Centrold Calculation Defuzzificatlon Method 



Another method for calculating the output value is the centroid method. In this method, a weighted average 
of all the active rules determines an output by summing all of the applicable output variables over their 
relative membership values. Although this method is more computationally intensive, it creates a distinct 
output value based on the relative memberships of all of the active rules that apply (see Figure S). This 
method eliminates the prd>lem of multiple solutions observed with the maximum method. A processor 
architecture with a hardware multiply-accumulate featiffe like that of the TMS320 DSP family excels at 
this method. 

Figure S. Fan_Speed Output Centroid Calculation 
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Summary 

By using fuzzy logic, you can simplify complex control problems that once required a high-powered micro- 
processor to execute in real time; you can now execute them on a low-cost Texas Instruments TMS320 DSP 
or TMS370 microprocessor. The following application note shows the benefits of omtrolling a simple DC 
motor with fuzzy logic using a TMS320C14 digital signal processor. 



Implementation of 
Fuzzy Logic Servo Motor Control 

on a 
Programmable 
Texas Instruments TMS320C14 DSP 



Mathew George, Jr. 
Digital Signal Processing — Semiconductor Group 
Texas Imtrummta Incorporated 



Abstract 



This paper describes the implementation of a fuzzy logic compensator on a Texas Instruments TMS320C14 
DSP-based servo motor control development system. The system contains a real motor that is controlled 
by the programmable DSP. An on-chip debugger and servo motor program allowed both simple code modi- 
fication and interactive control of the motor. A fuzzy logic algorithm was directly substituted for the origi- 
nal PID algorithm; this resulted in comparable motor response and algorithm performance. This imple- 
mentation proves ^e feasiliify §imsiMtm tasEZ^ lQ0:c-bm^<0Btm mmx &mM m a real system. 



Introduction 

Fuzzy logic is relatively new theory. Most of the readily available hands-on fuzzy logic system examples 
have been software simulations or bulky real systems. However, a TI commercial microprocessor (the 
TMS320C14) can serve as a simple, real-time, real-system platform for applying and investigating fuzzy 
logic. This facilitates both the understanding and implementation of fuzzy logic as a real-time program- 
mable solution for the general engineering public. 

Servo motor control is a viable and useful implementation. A programmable PID motor control board uses 
a Texas Instruments TMS320C14 chip and has an actual motor whose performance can be observed. Faster 
and newer parts are available, but the 'C14 is optimized for motor control with such features as on-board 
pulse-width modulation (PWM) generation capabilities. The board also has on-chip debugger code and 
interactive PID control code. The PID compensator code (written in TMS320 assembly language, which 
is upward-compatible with code executed by such newer fixed-point TI DSPs as the TMS320C25 or 
TMS320C50) uses position and velocity of the motor for inputs and motor input current as the output. The 
control code will allow the operator to interactively change such values as servo position and velocity and 
to monitor position enor. Modifying the code for fiizzy logic reguiredi^plaGing the PID compensator sec- 
tion of the code with a fuzzy logic compensator. 

The membership function for the compensator defines the error between present motor position and the 
desired (command) position of the controller. Five linguistics variables characterize the function: negative 
medium, negative small, zero, positive small, and positive medium. TTie function is fepresented as overlap- 
ping isosceles right triangles for ease of fuzzification. Various rules for an inverted pendulum control sys- 
tem (ball and stick) are explained in [1]. These same eleven rules were used for servo motor control. 

The algorithm's three sections are implemented as a series of software loops: fuzzification (i.e., input evalu- 
ation), fuzzy inference (rule contribution that uses a table look-up), and defuzzification (using center-of- 
gravity method). It is based on an algorithm developed for [4]. Each section uses various arrays that are 
modified in that section. The sizes of these loops are directly proportional to the number of inputs, outputs, 
and rules used in the system. 



Servo Motor System (Power-14) 

The heart of the Power-14 board is a Texas Instruments TMS320P14 chip (one-time programmable 
TMS320C14). (See Figure 1). 
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. P^ira t. Power 14 SyttMfi 




You communicate to the board through an RS-232 serial port connection by using standard terminal or ter- 
minal emulation software (such as Procomm). The 'P14 peripherals are optimized for control applications. 
An event manager can be operated in a PWM (pulse-widdi modulation) mode that is ideal for motor control. 
Monitor code is loaded into 'P14 external memory and run, which provides a command line debugger. The 
debugger has all standard debugger functions, such as memory dumps and modification, stepping through 
code, breakpoints, etc. The monitor can be used to load and run the servo motor program with PID compen- 
sator. The program is interactive and lets you control the motor from the keyboard. Position, velocity, and 
the PID values can be set. Data acquisition functions allow an ASCII text input stimulus table to be loaded 
onto the Power- 14, an acquisiticm run to be executed, and the resulting ASCII output table to be written 
to a PC file for graphing. The rest of the board contains support circuitry: amplifiers for the motor and a 
serial port interface. An encoder on the motor is used as a position sensor for a compensator input. The 
velocity is found by executing a back-diffeience. Note that there is no separate velocity sensor. 



PID Imptomefitation 

The source code for this system implements the PID compensator in one file (See Appendix A). The algo- 
rithm is a direct implementation of the PID equation. The Proportional, Integral, and Differential variables 
are derived from the motor encoder sensor detectiAg position. On each cycle, the present position is taken 
from the encoder and stored in Position. The error is found by subtracting Position fi:om DesiredPosition 
and storing it in ErrNow. Thus, ErrNow is the position input for the proportional section of the PID. A back- 
difference is then taken with ErrNow and ErrLast (the error in the previous cycle) to give ErrDiff . ErrDiff 
is an approximation of the velocity and is tfaerdbre die second input (Differential) of the PID. The Integral 
is found by adding the EnNow value and storing it in Kintegrator. The following equation tiien balds the 
compensator output: 

PWM = 9* P* ErrNow + 9 * I * Kintegrator + 9* D* ErrDiff (i) 

The value is scaled for the PWM mode and stored in NewServo. Thus, for the compensator output, the 
actual motor current input value is converted to a PWM value for implementation. The PWM output is then 
sent to power amplifiers that finally drive the motiH'. The PWM frequency can be controlled firom the PID 
program. 
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Fusy Lo^thmiry for Servo M^r Control 



The membership function for this system is simple. Five fiizzy lo^ ranges (linguistic variables) were cho- 
sen to remain consistent with [1] and used for both inputs. Sets of values are represented as overlappii^ 
isosceles right triangles (Figure 2). 

Figure 2. Membership Function 



















/ 




V 


4 




V 




\ 


-3 


-2 


-1 





1 


2 


3 



Key 

PMspositive medium 
PSsposHive small 
Z-zeio 

NS>negative smaU 
NM>negative medium 

TWo input variables, position (Theta) and velocity (dTheta) of the motor, are used in this fuzzy logic system. 
One output variable, motor current, is used, which will be proportional to the PWM output that is actually 
written. The rules for the compensator as mentioned were taken from [1]. Thus Theta, dTheta, and motor 
current operate according to the following "if a and b, then c" rules, as shown in Table 1. 

Table 1. List of Rules 



If Theta s 


And dTheta s 


Then Motor Current s 
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These rules can then be indexed according to the sesh ^bowa^ Figure 2. The mapping seen in Ikble 2 
will be used in TMS320 programming. 

Table 2. Indexed List of Rules 



If Theta = 


And dTheta = 


Then Motor Current = 











1 





-1 


2 





-2 


-1 





1 


-2 





2 





-1 


1 





-2 


2 
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-1 
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-2 
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-1 





-1 


1 






The defiizzification is done by the center-of-gravity method. 

Fuzzy Logic Implementation 

Arrays are used in the fuzzy logic calculations and modified in the various loops that implement the com- 
pensator. Appendix B lists the TMS320C14 code. Figure 3 shows the three sections of the compensator: 
fuzzification, fuzzy inference, and defuzzification. Notation for the aaays follows C language standard, 
with subscripts from to n-1. The 'C14 algorithm is based on an algorithm developed for [4]. The figure 
key summarizes the values of the system that will be used in the examples in tliis section. 
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RgtwieS. listed Arrays 
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for i=0 to ninputs-1 
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X[nlnputs][nmembs] 
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Fuzzy 
inference 


for i=0 to nrules-l 
for j=0 to ninputs-1 
for I<bO to noutputs-1 








U[noutputs][nmembs] 


Defus^cation 


for i=0 to nmembs-1 
for i=0 to nmembs-1 






u [outputs] 

r 





FU2zy Logte implementatlort as a Series of L^ps 



toy 
nmembs=5 
ninputs^ 
nrulessll 
noutputs^l 

I, j, l< = loops within loops 

Note which arrays are modified in each section and their size boundaries. They are discussed later in more 
detail. The two inputs seen in x[ninputs] are position (found from the encoder) and velocity (found from 
an approximation of the derivative by taking the back-difference of the position). These are the Theta and 
dTheta variables described previously. As the compensator code begins, the position and velocity inputs 
(ErrNow and ErrDiff, respectively) arc copied to the array X[ninputs]. The position is mapped so that one 
rotation of the motor ranges from -255 to +255 (See Figure 4 a). This relationship is mapped onto the x-axis 
of the membership function and thus fiizzifies the position of the motor (Figure 4 b). Note that this figure 
is not drawn to scale. 
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Figure 4. Motor Shaft/M«mb«3htp Function iUlapfiIng 
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In the fuzzification loop, the degree of membership of each input relative to the input membership function 
is evaluated and written in the array X[ninputs][nmembs]. The value for a particular linguistic variable is 
the y value of the triangle for a particular value of x (See Figure 5). 

Figure 5. Analysis for One Linguistic Variable 




15 



The y value is found bjr iidng the sinqde algebraic equation for a line (y=mx+&). This equation may be 
geometrically reduced to one of die two foUowing equations, depending on whidi side of the triangle the 
value of X lies: 

if {x :xO < X < x\} then y = (x - xO) / (xl - xO) else (2) 
if{x:xl<x< x2} then y = Qc2 - x) / (x2 - xl) 

The division needed is costly on most microprocessors, usually requiring at least a number of cycles equiv- 
alent to the number of bits of the number being divided. But if the slope m is made equal to 1 by causing 
the elements of the membership functions to be isosceles right triangles, the equation can be reduced to 
y=x-xO. This translates into a simple one-cycle subtraction. Note that forcing the elements of the member- 
ship functions to be isosceles right triangles also forces the peak of the membership function to no longer 
be 1. Rather the peak value^i-xOsx2-xisgS (as seen in Example 1). This action also eliminates the need 
for ming a Q format [6] to represent the fractional values from Equation 2 if the triangle were not isosceles. 



Example 1 demonstrates this fuzzification calculation. If the position input value were 100 (i,e., x[0]=lDO), 
it would have nonzero degrees of membership in the PS and PM linguistic variables. You can also see this 
in Figure 5. To calculate the actual degree of membership value, the value for;c is plugged into Equation 
2 for both PS and PM boundaries. Thus, in PS the contribution is 70, while in PM it is 15. The rest of the 
linguistic variables are because there is no contribution. (Of course, in software, all linguistic variables 
must be evaluated.) For this example, X[0][nmembs]=[0, 0, 0, 70, IS]. 

Exampte 1. Fuzziflcation Example 




-255 



x=^00, thus 
y=70andy =15 
PS PM 



The next step involves fuzzy inference. In this loop, the maximum and minimum functions, as explained 
in [2], are implemented. In the actual code, the array indexing of the membership values is made nonnega- 
tive by adding a bias of three. Therefore, instead of NM to PM being indexed from -2 to +2, as seen in 
Figure 4, they are indexed &om to 4. Table 3 shows how the 11 rules are indexed and reindexed in the 
array RULE_TABLE[rule][input-K»i^t]. 



16 



Table 3. Original and Relndexed Table of Rules 



Original Index 


Reindexed 


x(0) 


x(l) 


u(0) 


x(0) 




u(0) 











2 


2 


2 


1 





-1 


3 


2 


1 


2 





-2 
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-1 
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1 


2 
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-2 
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-1 


1 


2 


1 


3 





-2 


2 


2 
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1 


-1 


2 


3 


1 





1 


-2 


2 


4 





1 


-1 





3 


1 


2 


-1 


1 





1 


3 


2 



Some explanation is required for this decoding. The indexed rules match the explicit rules as described in 
Table 1. Also, the values given by accessing the RULE_TABLE are limited to the values indexed by 
nmembs. This characteristic is heavily used in the index manipulation and allows nmemb to be inter- 
changed with RULE_TABLE[rule][input-K>utput] in the appropriate parts of the algorithoi. 

To find the minimum value of X[ninputs] [nmembs] decoded from the rule table inputs and stored in minZ 
(which is initialized to the maximum y value — in this case, 85), the equation is: 

miwZ = min (X[input] [RULEJTABLE [rule] [ityjut]]) p) 

(Since noutputs=l, the loop is simplified, and minZ does not need to be the general case array 
minZ[nrules]). Note that only the first two columns (the input columns) of RULE_TABLE are used in this 
part of the fiizzy inference section. 

Then, for each rule, the max is taken of the output value U[nmembs], which is initialized to 0. The general 
case U[output][ nmembs] is simplified because only one output is decoded from the rule table outputs and 
minZ. This equation can be summarized as: 



U[nmembs] = U{RULE_TABLE[rule] [ninput + output]) 

= max (U[RULE_TABLE[rule] [ninput + output]], minZ) 




Note that in this section only the last column (the output column) is used. The following equation summa- 
rizes the minimum and maximum functions that are executed for each rule to result in the array U[nmembs] 
by plugging Equation 3 into Equation 4: 

U[nmafU>s] = max {U[RULE_TABLE[rule] [ninput + output]] , 
min Q[[input] [RULE_TABLE[rule] [ir^ut]])} 
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The following example iliOSteaMi v^m&s ieation. Ofie cjKde of the loop for rulesO is shown. 

For the input array: 



X[ninputs] [nm&nbs] = I 70 15 | 
= 1 60 25 I 



(6) 



and for rule=0, input=0, and inpat=l, plug into Equation 3: 

minZ = min \ {X[0] [2]) | = min \ 70 \ =25 
I iX[l] [2]) I I 25 I 



(7) 



then for the max nile=0, inputsQ, and input=l, plug into Equation 4: 

U[nmembs] = U[2] = max \ U[2] | = max \ | = 25 

I minZ I I 25 I 



(8) 



This process continues for the list of 11 rules so that the array U[iunembs] contains the maximum values 
of the min — ^i.e., the contribution of each rule to the inference. 

The final step involves defuzzifying the U[noutputs][nmembs] array. Since noutputs=l, U simplifies to 
U[imiembs]. The U[nmembs] array now has five values in it for its corresponding five positions. The cen- 
ter-of-gravity calculation is done with two loops. The first finds the numerator by using the multiplier to 
weight U[nmembs] by its position. The second loop finds the denominator by summing the position. The 
inevitable 16-bit divide loop then finds the output value u[output]. This u[output] represents the motor cur- 
rent mentioned mFuzzy Logic Theory for Servo Motor Control (page 12). The divide operation is done on 
the 'C14 by a 16-cycle loop. This value is then scaled for the output and written to the NewServo memory 
location ^t sends it to ihs PWM j^nerator. 



The fuzzy logic code with two inputs, eleven rules, one output, and five linguistic variables requires about 
2000 instruction cycles to execute on the TMS320C14. This means a 400-^8 period (2.5-kHz frequency) 
because one instruction cycle is 200 nanoseconds on the 'C14. The update period of the motor is 3.8 ms 
(263-Hz frequency), so the fuzzy code is obviously adequate (for comparison, the PID code execution 
required nine us (111-kHz frequency), for the update period). Figure 6 shows the PID performance for a 
step with the classical control overshoot and settling times. 



As an example of defteification, 
if U[mnemb] = [0 15 70 35 OJ, 
then the output is: 

U[output] = 0*(-170) + 15*(-85) + 70"^{0) + 35*(85) + 0*(170) 



= 14.17 



These loops thus evaluate the control input necessary for the servo motor. 



Results 
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The update period is adequate for applications such as servo motor, robotics, motion control, and automo- 
tive control. Better performance may be needed for such applications as hard-disk drives. Later-generation 
digital signal processors, such as the TMS320CSx, operate at up to 25 nanoseconds with much mem e£B- 
cient instructicm code. 



Conclusion 

The final fuzzy logic system behaved favorably when compared to the conventional PID system. The sys- 
tem proved the feasibility of implementing a real-time fuzzy logic servo motor control. Proof of the fuzzy 
control was shown by positioning the motor spindle with an error outside the membership function, thus 
causing the control to desist. Further development could include a graphics display and the ability to vary 
the rules. The TMS320C14 board fits in a 12 x 8 x 6-inch suitcase conveniently and requires only an AC 
power supply and an RS-232 keyboard coaaecti(»i. This product easily demonstrates real-system fuzzy 
logic control on a microi»ooessor. 

Fuzzy logic has great potential as a programmable solution for general engineering. For applications where 
performance is a priority, a hard-wired silicon solution (which may even configure as a microprocessor 
peripheral) based on [4] is being devel(q)ed. The programmable solution may assist in ttie transition to this 
hard-wiied opticMi, dqtrading cm the sofltware/hardwaie tradeoffs. 
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APPENDIX A 
PID Code 



title "Servo compensator $Revision: 3.4 $" 

*••••••••*••••*••*•*«********•****••*•***•***••**************••<***••••*••*• 

SHeaden: C:/src/cl4/ps/vc5/comp^_v 3.4 01 Oct 1991 17:24:18 "$ 

$config$="A'8 /K! /L;*j«f.def /R:*-/B8(r 
!coiifig!="/Mcomp.s'' 

INAMEI 

comp.s 

'EAIHS! 

modules 

!0! 

DESOUPnON 

Servo compensator - uses a PID algorithm 

PRINCIPLE AUTHORS: 
Dave Sewhuk 

CREAnONDATE: 

December 23, 1990 22:54:08 

COPYRIGHT NOTICE: 

(C)Copyright 1990 Teknic Inc. All rights reserved. 

!end! 

HEADERS imUZED 

.include "macrodef.inc" 
.include "cl4io.inc" 

******** **«««a««*«««********4>*«**«*««*4iiii***««*******»***«iit«***«iii****iii«iii*«** 

INAME! 

comp.s 



IPATHS! 
•0! 



Exported Variables 

.def PwmChannel^PwmPeriod 
.def NewServo3x^DifF,Kintegrator 
.def EirNow3rrLastJ)esiiedPo^c»i 



.bss PwmChannel,! 

.bss PwmPeriod,! 

.bss NewServo.l 

.bss ErrDiff,! 

.bss Kintegrator,! 

.bss ErrNow.l 

.bss ErrLast,! 

.bss DesiredPositioii,! 



; Current PWM channel 

; Current period of PWM channels 

; New servo value 

; Error difference 

; Error integrator 

; Current error function 

; Last error function result 

; Servo to this position 



!end! 



INANfE! 

comp.s 

IPATHSI 

Local Variables 

!0I 
!end! 



.text 

.def servo_iSiR 

; 'NAME! 
servo_ISR: 

• IRATHS! 



This is the servo compensator. The distribution algorithm is 
the PID format. The compensator performs the following ftinctiiffia 
PWM = 9*P*error + SUM(9*ek)*I + 9*(err difO*D 
with the PWM output clipped to timer hardware limits. 



No changes to foreground environment. 



STACK LEVELS: 


EXTERNAL DATA REFERENCES 



.ref ONI^CONSTFFFF ; Useful c®mt^ts 



fimctionsNall 

ftmctions\intemipt 

modules\utils.s 



roi 

DESCRIPTION: 



RETURNS: 



.ref Position 



; Position Irom feedback ISR 
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.ref Pgain,Igain,Dgain ; Servo constants fircm user interface 
.ref ISRJTMP ; Scratch 



!skip start! 



; Saturate math 



SOVM 

1 

LAC ErrNow 
SACL 



ZALH 
; Calc new enor 

SUBH 
SACH 



BrLast 

DesuedPosition 

Position 
^Now 



Run PID algorithm on new sample 



LAC ErrNow 
SUB ErrLast 
SACL 

ZALH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
SACH 



ZALH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
ADDH 
SACH 

; Saved scaled up error 

ZALH 
; Update integrator 

ADDH 
SACH 



ErrDiff 

ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 
ErrDiff 

ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ErrNow 
ISR TMP 



Kintegrator 

ErrNow 
Kintegrator 



; Save old eifor 



; Update error difference 



Scale up value a bit/saturation 

»2 

•3 

•4 

•5 

*6 

•7 

*8 

•9 



; Scale up/saturation 

;*2 
;*3 
;*4 
;*5 
;*6 
;*7 
;*8 
;*9 



LT Pgain 
MPY ISR_TMP 
PAC 

LT Igain 
MPYKintegrator 

LTA Dgain 
MPYErrDiff 
APAC 
SACH 



NewServo 



; Do proportional part 

; Do integral part 
; Do differential part 
; Save result 



Check the answer so that it fits withm the limits of the timer's 
period: [O..PwmPeriod*4) 



pwmisr20: 



LAC PwmPeriod,! 
SUB 0NE,3 
SUB NewServo 
BGZ pwmisr20 
LAC PwmPeriod,! 
SUB 0NE,3 

SACL NewServo 
B pwmisrSet 

LAC NewServo 

ADD PwmPeriod,! 

BGEZ pwraioSet 

ZAC 

SUB PwmPeriod,! 

SACL NewServo 



Set Output PWM to new value 
pwmisrSet: 

LAC PwmPeriod,! 
ADD NewServo 
SACL NewServo 
.if ChipV!R! 
.ref C0NST3 
SUB C0NST3 
BLEZ pwmisrlO 
; No fixes needed 

LAC NewServo 
ANDC0NST3 
SUB C0NST3 
BNZ pwmisr!0 
LACK 4 
ADD NewServo 



Oieck positive: 2*Period-4 
- value 

Oops, pin at biggest 

; Ne^tive OK? 
; Add largest offset 

; Make largest negative 



; Get center period value 
; Add newly computed value 

; PWM bug workaround 



; Lower 2 bits set? 
; OK, No fixes needed 
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SACL 



NewServo 



pwmisrlO: 



.endif 

LACK ActionBank 

SACL ISR_TMP 

OUT ISR_TMP,BSR 

OUT NewServo,ACTO 

LAC PwmPeriod,2 

SUB 0NE,2 

SUB NewServo 

SACL NewServo 

OUT NewServo^CTl 

RET 



; Iskip end! 
; !END! 



;ENDOFFILE 



; Set bank for timers 



; Output channel 1 timer 

; Calculate complimentary output 



; Set complimentary chaimel 2 timer 
; All done., back to Hkm ISR already 
; in progress. 



.end 
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APPENDIX B 
Fuzzy Logic Code 

.title "Fuzzy Servo compensator Revision: 3.4 $" 

SHeader:: C:/src/cl4/ps/vcs/comp.s_v 3.4 01 Oct 1991 17:24:18 "$ 

$config$="/T8 /K! /L;*.ref.def /R:*-/B^" . 
iconfig! ="/Mcomp.s" 

!NAME! 
compfuz.s 

lEATHS! 

modules 

!0! 

DESCRIPTION 

Servo compensator - uses a fuzzy algorithm (modifying PID code) 

PRINCIPALAUTHORS: 
Dave Sewhuk 

Joe George w/ Fuzzy based on C++ code by Phillip Thrift 
CREATION DATE: 

December 23, 1990 22:54:08 

Jan. 1992 
COPYRIGHT NOTICE: 

(C) Copyright 1990 Tbkmc Inc All rights reserved. 

(C^ Cof^gbt 1992 Texas Instruments Incwporated. All rigjits r^rved. 

!end! 

HEADERS mUJZED 

.include "macrodef.inc" 
.include "cl4io.inc" 

INAME! 

comp.s 



IPAraS! 
!0! 



Exported Variables 

.def PwmChannel,PwmPeriod 
.def NewServo,ErrDi£f,Kinte^ator 
.def &rNow,EfrLast,DesiredPosition 
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.def ninput, noutput, triangle, nmemb 

.def NM, NS, ZE, PS, PM, nloc, NMw, NSw, Zw, PSw, PMw 

.def nrule, rule_table, input, output, rule, memb 

.def l«ic,mim, den, TEMP, TEMPH,x,X,U 



.bss PwmChannel.l 

•bss PwmPeriod,! 

.bss NewServo,! 

.bss ErrDiff,! 

.bss Kintegrator,! 

.bss ErrNow.l 

.bss ErrLast,! 

.bss De^dB^osidon,! 



; Current PWM channel 

; Current period of PWM c^nnels 

; New servo value 

; Error difference 

; Error integrator 

; Current error function 

; Last error function result 

; Servo to Ms position 



; The following are fuzzy section, initialized and uninitialized 



scalein .set 256/3 

scaleout .set 2*100/3 

; .sect "pfuzzin" 

.asect "pfuzzin", lOh 
.label prule_table 



rule table: 



.word 2,2,2 
.word 3,2,1 
.word4,2,0 
.word 1,2,3 
.word 0,2,4 
.word2,l,3 
.word 2,0,4 
.word 2,3,1 
.word 
.word3,l,2 
.w<»dl3>2 



2,4,0 



ninput 
noutput 



.word2 
.wotdl 



triangle .word 3 

nmemb .word 5 

NMaddr .word NM 

NM .word -3*scalein, -2*scalein, -l*scalein 

NS .word -2*scalein, -l*scalein, 0*scalein 

2E .word -l*scalein, 0*scalein, l*scalein 

PS .word 0*scalein, l*scalein, 2*scalein 

PM .word l*scalein, 2*scalein, 3*s(»lein 



nloc 

NMwaddr 

NMw 

NSw 

Zw 

PSw 

PMw 



.word 5 
.word NMw 
.word -2*scaleout 
.word -l*scaleout 
.word O*scaleout 
.word I'scaleout 
.word 2*scalequt 



nrule 



.word 11 



rule tableaddr .word rule table 



input 

output 

rule 

memb 

loc 

num 

den 

minZ 



.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 



TEMP 

TEMPH 

0NE2 



.word 
.w<xd 
.w<»d 1 



xaddr 

Xaddr 
Uaddr 



.word X 

.word X 
.word U 



rule_end: 
rambeg 

X 

X 
U 



.used "fuzdn*', nile_end-^e_table 

.usect "fuzzunin", 2 
.usect "fuzzunin", 2*11 
.usect "fuzzunin", S 



; !end! 



!NAME! 



lEATHS! 



comp.s 



Local Variables 



!0! 



; lend! 
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.text 

.def servo ISR 



; !NAME! 
servo_ISR: 

IPATHS! 



functions\all 

ftmctio]lS^intert^ 

modules\utils.s 

!0! 

DESCaUPnON: 

This is the servo compensator. The distribution algorithm is 
mocked fma PH) to fuzzy format. 

RETURNS: 

No diaoges to foreground environroent. 

STACK LEVELS: 


EXTERNAL DATA REFERENCES 

.ref Position ; Position from feedback ISR 

•ref ONE,CONSTFFFF ; Useful constants 

.ref ISR_TMP ; Scratch 

iskip start! 

SOVM ; Saturate math 

; First need to copy data over from program 

LACK prulejable 
ADD ONE, 11 
ADD ONE, 8 

LDPK 
LARK ARO, 6Fh 

LARK ARl, rambeg ; beginning of RAM block 

RAMLOOP:MAR M 

TBLR •+ 
LDPK 1 
ADD ONE 
LDPK 
MAR»,0 

BANZ RAMLOOP 

; Get Heta and dUieta valties 



LDPK 

LAC ErrNow 
SACL 

ZALH 
SUBH 
SACK 

LDPK 

MAR 
LAR ARO^ddr 
SACK *+ 

LDPK 
LAC ErrNow 

SUB ErrLast 
SACL ErrDiff 

LDPK 
SACL 



ErrLast 

DesiredPosition 

Position 

ErrNow 



•.0 







; Calculate X 

; indices for X[input][memberfii] 
x&: 

ZAC 
SACL 
xfoil: 

SACL memb ;forQsO; 
LAC 
SUB 
BLEZ 
xfriml: 
SUB 
BLEZ 
LAC 
ADD 
SACL 
LAR 
LAC 
LT 
MPY 
APAC 
SACL 
LAR 
CALL 



SACL 

LT 

MPY 

PAC 

ADD 



input 
ZAC 
memb 
ninput 
input 
ilend 

LAC mnemb 
memb 
mlend 
xaddr 
input 
TEMP 
ARO.TEMP 
NMaddr 
triangle 
memb 

TEMP 
ARl.TEMP 
membaship 

TEMPH 

nmemb 

inptit 



; Save old error 
; Calc new error 

; store x input value (Theta) 



jdTheta 



; for (i=0; 

; i<mnput 

; j<nmemb 

; load &x[input] ARO 

; load &memb into ARl 



; Get membership value with AR's loaded 
; and leave answer in ACC 
; Save value 

; Calculate index for X[iiiput}[memb] 
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ADD 


Xadcb 


SACL 


TEMP 


MAR 


M 


LAR 


ARl.TEMP 


LAC 


TEMPH 


SACL 


• 


LAC 


memb 


ADD 


0NE2 


SACL 


memb 


B 


xfmnl 


mlend: 


LAC input 


ADD 


0NE2 


SACL 


ini»it 


B 


xfoil 


ilend: 


NOP 


; exit input loop 


; Zito out U values 


LAR 


ARO, nmemb 


LAR 


ARl, Uaddr 


ZAC 




Uzl:MAR 


M 


SACL 


*+ 


MAR 


•.0 


BANZ 


Uzl 


; calculate value of Z 


zSa: 




ZAC 




SACL 


rule 


zfiirl: 


ZAC 


SACL 


input 


LAC 


nrule 


SUB 


rule 


BLEZ 


zrlend 


; LAC 


ONE?, 


LAC 


PM 


SACL 


minZ 


zfhil: 


LAC ninput 


SUB 


input 


BLEZ 


zilend 


LAC 


ninput 


ADD 


noutput 


SACL 


TEMP 


LT 


TEMP 


MPY 


rale 


PAC 




ADD 


input 



; Store value in X[input][memb] 



;m«iib+4- 



; input++ 

; exit memb loop 



;for(i=0; 

;forG=0; 
; i<nrule 



; initialize Z[rule] 

; kluge since PM = l*scalein 

; i<niaput 



; find &RULE_TABLE^6][input] 



ADD 


rule_tableaddr 


SACL 


TEMP 


MAR 




LAR 


ARO.TEMP 


LAC 




SACL 


TEMP 


LT 


nmemb 


MPY 


input 


PAC 




ADD 


TEMP 


ADD 


Xaddr 


SACL 


TEMP 


LAR 


ARO,TE»ilP 


MAR 


*.o 


T AC 


ininZ 


SUB 


* 


BLZ 


notmin 


LAC 


• 


SACL 


minZ 


notmin: 


LAC input 


ADD 


0NE2 


SACL 


input 


B 


zfiiil 


; chedc max for U 


zOend: LAC 


ninput 


ADD 


noutput 


SACL 


TEMP 


LT 


TEMP 


MPY 


rule 


PAC 




ADD 


ninput 


ADD 


rule_tableaddr 


SACL 


TEMP 


MAR 


*.o 


LAR 


ARO.TEMP 


LAC 


* 


ADD 


Uaddr 


SACL 


TEMP 


MAR 


*.o 


LAR 


ARO,TEMP 


LAC 


minZ 


SUB 


• 


BLZ 


notmax 


LAC 


minZ 


SACL 


• 


notmax: 


LAC rule 


ADD 


ONE2 



; find RULE_TABLE[rule][input] 



; TEMP = RULE TABLE[rule] [input] 
;find&X 



; &(X[input][RUl£_TABLE[rule][input]) 
; check for min 



; input++ 



; get &U[RULE_TABLE[rule][input+output] 



; total kluge for input+output 



; got RULE_TABLE[rule][input+output] 



;U[RULE_TABLE[rule][input+output]] 
;if U<minZ (looking for max) 

store new U 

; rule-H- 



as 



SACL rule 
B zfnrl 
zrlend: NOP 



; Need to defuzzify output, u = num/den = sigma U*Loc/ sigma U 
; (ignoring wt.) 
defuzz: 
ZAC 

SACL loc ; for (1=0; 

SACL 
SACL 
LAR 
LAR 



; Sigma loc*U 

dfll: 
SUB 
BLEZ 
MAR 
LT 
MAR 
MPY 
ZALH 
ADDS 
APAC 
SACK 
SACL 
LAC 
ADD 
SACL 
B 

dlendl: LAC 
SACL 
NOP 

; Sigma U 

ZAC 

SACL 

SACL 

LAR 

MAR 
dfl2: 

SUB 

BLZ 

LAC 

MAR 

ADD 

SACL 

LAC 



loc 

TEMP 
TEMPH 
ARO, NMwaddr 
ARl, Uaddr 

LAC nmemb 
loc 

dlendl 

*,0 

*+ 

M 
*+ 

TEMPH 
TEMP 

TEMPH 

TEMP 

loc 

0NE2 

loc 

dfll 

TEMP 

num ; Store 



loc 

TEMP 

ARO.UaAlr 

*,0 

LAC nloc 
loc 

dlend2 
TEMP 
*,0 
•+ 

TEMP 
loc ; ++i 



; i < nloc; 



; for (1=0; 



; i<nloc; 



ADD 


ONE2 


SACX 


loc 


B 


dfl2 


dlend2: LAC 


TEMP 


BZ 


dzero 


SACL 


den 


NOP 




B 


cdiv 


dzero: 


ZAC 


B 


endfuz 



; divide num/den and scale for new seivo 



cdiv: 



LT 

MPY 

PAC 

SACK 

LAC 

ABS 

SACL 

ZALH 

ADDS 

ABS 

LARK 



DIV: 



BANZ 

SACL 
LAC 
BGEZ 
ZAC 

SUB 
SACL 

done: 
LAC 
B 

membership: 
MAR 
LAC 
MAR 
SUB 
BLEZ 
MAR 
LAC 
MAR 
SUB 
BLZ 
MAR 



MAR*,0 
TEMPH 
den 

TEMP 
den 

den 

TEMPH 
niim 

AR0,1S 

SUBC 
DIV 

TEMPH 

TEMP 

done 

TEMPH 
TEMPH 

ZAC 

TEMPH 

endfuz 



*,0 
* 



M 



'+ 

zvalue 

•.0 
* 

M 
•+ 
divl 
♦.0 



;Sign 



den 



; Quotient 
;Sign 



; New Servo 



; &x[input] 

; &memb 
; if x>xO 



;if x>xl 
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LAC 


• 




MAR 


M 




SUB 


•+ 


;ifx>x2 


BLZ 


div2 




B 


zvalue 




divl: MAR 


•,ARO 


;x-xO 


LAC 


* 




MAR 


•,AR1 




MAR 






MAR 






SUB 


* 




RET 






div2: MAR 


•,AR1 


;t2r-x 


MAR 


*_ 




LAC 


• 




MAR 


•,ARO 




SUB 


• 




RET 






zvalue: LACK 





; value = 


RET 






endfuz: LDPK 


1 




SACL 


NewServo 


; Save result 


ZAC 




; negatioa Istx^ 


SUB 


NewServo 




SACL 


NewServo 




; yoda: B 


servo_ISR 





Check the answer so that it fits within the limits of the timer's 
period: [O..PwmPeriod*4) 



LACK 
SACL 
LAC 
SUB 
SUB 
BGZ 
LAC 
SUB 
SACL 
B 

|iwmisr20: 
LAC 
ADD 
BGEZ 
ZAC 



OFFh 

PwmPeriod 

PwmPeriod,! 

0NE,3 

NewServo 

pwmisrZO 

PwmPeriod,! 

0NE,3 

NewServo 

pwmisrSet 

NewServo 

PwmPeriod,! 

pwmisrSet 



Check positive: 2*Period-4 
- value 

Oops, pin at biggest 



; Negative OK? 
; Add largest ofiiset 



SUB PwmPeiiod,! 
SACL NewServo 



; Make largest ne^tive 



; Set Output PWM to new value 
pwmisrSet: 



LAC 


PwmPeriod,! 


ADD 


NewServo 


SACL 


NewServo 


.if 


ChipVlRl 


.ref 


C0NST3 


SUB 


C0NST3 


BLEZ 


pwmisrlO 


LAC 


NewServo 


AND 


C0NST3 


SUB 


C0NST3 


BNZ 


pwmisrlO 


LACK 


4 


ADD 


NewServo 


SACL 


NewServo 


pwmisrlO: 




.endif 




LACK 


ActionBank 


SACL 


ISR_TMP 


OUT 


ISR^TMP.BSR 


OUT 


NewServoACTO 


LAC 


PwmPeriod,2 


SUB 


0NE,2 


SUB 


NewServo 


SACL 


NewServo 


OUT 


NewServo^CTl 


RET 





; iskip end! 
; !END! 



;ENDOFnLE 
.end 



; Get center period value 
; Add newly computed wato 

; PWM bug workanxmd 



; No fixes needed 
; liower 2 bits set? 



; OK, No fixes needed 



; Set bank for timers 



; Output channel 1 timer 

; Calculate oimplimentary output 



; Set complimentary channel 2 timer 
; All done., back to the ISR already 
; in progress. 
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The Programmable 
Fuzzy Logic Array 



Philip Thrift 
Central Research Laboratories 
Texas^Mx^enis Incoiporafed 



introduction 



The Programmable Fuzzy Logic Array (PFLA) produces nonlinear multidimensional mappings by encod- 
ing fuzzy rule systems in a programmable array architecture. It can be used as a component of a graphical- 
user interface (GUI) for designing fiizzy controllers, as well as a software blueprint for mapping onto VLSI 
hardware. An advantage of the PFLA over other fuzzy representations, such as &e FAM (Fuzzy Associa- 
tive Memory) [2], is the PFLA's ability to ea^y visualize several in|Hit8 and ou^ts simultaneously. 



Data Flow 



Figure 1 shows the PFLA data flow. Succeeding text describes PFLA in general mapping terms from P 
inputs to Q outputs, but only P = 2 and 2 = 1 are shown in Figure 1. Each item number in the text corre- 
sponds to a step of Figure 1. 



□ 
S 
H 



F^urel. PFLA Data Fi^ 

X2 



•x xkxx x" ^o ooo br 



I 



yy 



y 
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uennea over eacn inpui range [a;, Dij : i = i,.., are luzzy mcmoersnip luncuons f ... p"i, 

where tii is the number of membership functions defined for input {. Each membership function 
varies between and 1. The cases shown in Figure 1 are trapezoidal membership functions, 
which are defined in Appendix A. Other parametric families of membership functions can be 
substituted. Also, for each fuzzy membership function pi , there is a corresponding label. A typi- 

i 

cal labeling scheme (labels are not shown in Figure 1) for = 7 is: negative large (NL), negative 
medium (NM ), negative small (NS), zero (ZE), positive small (PS), positive medium (PM), pos- 
itive large (PL). Although the fuzzy sets for each input in Figure 1 appear symmetrically spaced, 
these are not necessarily the optimal settings. If trapezoidal memberships functions are used, 
four numbers = (i^^ (/^^ j must be specified for each membership function p. 

3. Each input xi is evaluated by each of its fuzzy membei^hip functions to produce a value y-/ : 

fi = F[ (x.) : i = 1,..., Pj = 1,..., n, (9) 

These values appear in the boxes as shown as a thermometer level, also refers to the box that 
contains its value. In the example shown, only two box^lor each input t«ive positive evaluation. 

4. Fuzzy rules are encoded in a oossbai pattern. (Corresponding to each of #ie input boxes in 

3, above, is a vertical wire dropping down. A horizontal wire crosses those vertical wires and 
also the vertical wires corresponding to output boxes in 6, below. A connection is indicated by 
a •. For each horizontal wire, there is, at most, one connection per input and output variable . Each 
horizontal connection pattern ^icodes a rule. For ^cample, the connection pattern on the first 
horizontal line encodes the rule: 

Ifxi is N$ aad Xj is ZE, Ui is PS. dO) 
Four rules are slio^ in Figure 1. ^ rules @an 1% specified by a table of numbers: 

ri , si ,...sj (12) 

where ^/ is in {1,...,«,,NULL}, is in {l,...,»i,-,NULL}. Here, mj is the number of fuzzy sets 
for output j. The kth horizontal wire specifies the rule 

F'i -^(fi G'^ (13) 

The NULL indicates that there is no cormec^n for this input/output variable. In Figure 1, the 
connections wot^d be 

3,4,5 
3,5,4 
4,6,3 
4,5,4 

5. Intercepting the horizontal wires between the inputs and outputs are the A ("wedge") boxes. 
A conventional A operator is the numerical minimum of the values (this is used in Figure 1), 
but other operators aie pc^sible (lor example, ia:oduct, or mf ei a set of so-odled t-aonns [1]). 
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At the hh A box, this is produced: 

h,= ^ (fl, /} fj) (14) 
If jj is NULL, this argument is omitted. 

6. The V (" vee") boxes are computed according to the values of the connections above them. A 
conventional V operator is the numerical maximum of the values (see Figure 1). The values of 
the V boxes are 

f j = ^ {hk : 4 = ;■] = l."-.j2 = (15) 

Other operators [1] can be substituted for this operator (for example, probabilistic sum: x®y 
= X + y-xy , etc.). 

7. Defined over each output range [c,-,4] : ' = h-,Q are fuzzy membership functions G\,..., Gj*' > 

where m,- is the number of membership functions defined for output i. Also defined for each out- 
put I is a vector 4° of locations quantizing the range: = ,„ /jij, where qi is the number of 

locatt^ sfM^iied for m^j^ i itf i@ttf 1, 17 to^^ti^ are sbmm fmi oa^iat 1. 

a ^ (l^, j = = 1,..., Mi s (16) 

Here, qj is applied by component to get the resulting vector. The ^ are precomputed and stored 
for each V box. Then, this is computed for output i: 

V, = V *(^» A • gi «^ A » ^ 

Here, the "wedge" and "vee" operators are not necessarily the ones used in the boxes above. 
In Figure 1, v/ is shown as a sequence of vertical bars; V * is the maximum operation, and A * 
is the product. 

8. The final stage is to compile i^- from vi for eadi output: 

II, = ^1-^ m 

Vi ■ 1 

where 1 = is a vector of Is. 

If, conesponding to eadi output fuzzy set ^ , ilmrs is a single distinct location /y^ = (j/^/^ ^ and A * 
is the pcodact operation, thet the above defiizzification procedure reduces to 

^ w] • g] • I] + ... + • g^' ■ r^' 

W] • g] + ... + H-;' • f^' 

This provides "weighted point mass c^tozificaticHi". In this case, only a "Wei^t" aad a location // 
must be stored for each V box. 

Figure 2 shows the state of the system for particular vapa\s xi,X2. Thermometer levels indicate the values 
in both the fuzzificaticm and the A and V boxes Qiere, min and max are the operations performed). 
Darkened rule connections indicate values that propagated through the array (since min and max are used, 
on each output line there is, in general, one max winner and one min winner). Figure 2 also shows the layout 
for a four-in^t, two-output (4-2) system. M aOUI, y%»u cm use point*and-cli(^IO Mst die rule connections 
and m^bershiptet^ioa posi^(»sig. 
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Figuroa^ PRJ^4>^ Layout 
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Summary 

This is the information that provides the setting for the PFLA: 

#Input fuzzy sets: ni,...,np 
#Output fuzzy sets: mi^..,mQ 
[ai,bi], .1 /"/,i = 1,...,P 

Operators: A V 



#Inputs: P 
#Outputs: Q 
Input fiizzy sets 

#Rules: R 
#Rule table: 

Output defiizziflers 



Appendix A 

A trapezoidal membership function on an interval [a,b] is specified by four numbers tx, t2, tj, 14, satisfying 



a s ti s t2 « tj £ t4 £ b 



The trapezoidal function is defined by 



TZ(a,ti,t2,t3,t4.bX: 





(x-ti)/(t2-ti) 
1 

(t4-xy(t4-t3) 




for X in [a,ti] 
for X in (ti,t2) 
for X in [t2,t3] 
for X in (t3,t4) 
for X in [t4,b 



Note that if tj = tz €>r ^ s (4, this part of the definition is void. 
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